Superfluid-insulator transitions of the Fermi gas 
with near-unitary interactions in a periodic potential 
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We consider spin-1/2 fermions of mass m with interactions near the unitary limit. In an applied 
periodic potential of amplitude V and period at,, and with a density of an even integer number 
of fermions per unit cell, there is a second-order quantum phase transition between superfluid and 
insulating ground states at a critical V = V c . We compute the universal ratio V c ma'l/h 2 at N = oo 
in a model with Sp(2A r ) spin symmetry. The insulator interpolates between a band insulator of 
fermions and a Mott insulator of fermion pairs. We discuss implications for recent experiments. 



An important milestone in the studies of ultracold 
atoms has been the observation of superfluidity in de- 
generate gases of fermionic atoms across a two-body 
Feshbach resonance [H, 12, 0, 0, @- As a function of 
the detuning, v, across the resonance, these systems in- 
terpolate from a Bose-Einstein condensate of diatomic 
molecules for large negative v, to a Bardeen-Cooper- 
Schrieffer (BCS) paired state of a Fermi liquid for large 
positive v. For v w 0, neither limiting description ap- 
plies, and we have a superfluid state of fermions with 
interactions near the unitarity limit. As has been em- 
phasized in recent work 0, 0, 0] this entire crossover has 
striking universal aspects, with all physical properties de- 
termined only by v and the density of the Fermi gas. 

In a separate development, ultracold gases of bosonic 
atoms were placed in an optical lattice potential 
With increasing lattice depth, the bosons exhibited a 
superfluid-to-insulator quantum phase transition. 

In the quest to realize strongly correlated quantum 
phases of interacting fermions, a recent experiment [101 ] 
has combined the techniques of the above experiment by 
studying fermions of mass m with near unitary interac- 
tions in the presence of an optical lattice potential of 
period ol and amplitude V. This paper shall demon- 
strate that the universality arguments can be extended 
to include the periodic potential after including a single 
energy scale associated with V. In particular, we find 
a universal phase diagram shown in Fig [1] below. This 
phase diagram has superfluid-insulator quantum phase 
transitions with a density of an even integer number (rip) 
of fermions per unit cell at critical amplitude V = V c 
which obeys 



a novel quantum state which is neither a band insulator of 
fermions, nor a Mott insulator of bosonic fermion pairs. 
Instead, it is in a interesting intermediate regime in which 
multiple single-particle bands are occupied. Previous 
computations of ultra-cold atoms in optical lattices have 

HQ, 
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relied on effective tight-binding models 
but such approaches are not expected to be quantita- 
tively accurate near the transition to the superfluid. A 
recent computation |15j | uses a method related to ours, 
but does not account for the off-diagonal couplings be- 
tween different reciprocal lattice vectors in Eq. (|6|) below; 
these are essential for a proper result, and make the com- 
putation much more demanding. Our 1/N expansion is 
also able to quantitatively account for the strong inter- 
actions between the fermions in the multiple bands. One 
consequence of the many occupied bands is that the com- 
putational requirements are demanding even at the lead- 
ing N = oo level. So our present numerical results will 
be limited to N = oo although we will set up a formalism 
that allows computations to all orders in 1/N. It is also 
worth noting that previous studies of ground state prop- 
erties [7[ found that 1 /N corrections were quite small in 
the unitarity limit, vwO. 

We consider 2N species of fermionic atoms ip^ (i = 
1 . , . N, a S {T,4}) coupled to a single field $ of s-wave 
Cooper pairs, or molecules. This is an Sp(2iV) gener- 
alization of the popular "two-channel" model, and the 
mysical case N = 1 can be accessed in 1/N expansions 
3,13]. The atoms experience an optical lattice potential 
V(r); we choose a simple cubic latice potential 



Vc = 



F nY (a h u) 



(1) 



where F nw is a universal function of a^v dependent only 
on the even integer rip. The transition at unitarity occurs 
a t F nF (0). We will show that F np can be determined in a 
1/N expansion in a model with Sp(2iV) spin symmetry. 
Explicit numerical results for the universal phase diagram 
and the function F nF will be presented below at N = oo. 

An important aspect of the physics accounted for by 
our analysis is that the insulator near the critical point is 



V(r) = V 



and our computations have simple generalization to other 
lattice structures. The density of the fermions is con- 
trolled by chemical potential /i and we assume that spin 
polarization is zero. The imaginary-time action of this 
many-body system, which includes all terms for a de- 
scription of the universal physics in the vicinity of a two- 
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body Feshbach resonance is: 
5 = / .lrd 3 
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Note that the Cooper pair field $ is actually a Hubbard- 
Stratonovich field that decouples the fermion interaction 
terms, and hence does not have a bare dispersion, or 
Berry's phase. Thanks to this, we have conveniently 
rescaled the $ field to absorb in it the interaction cou- 
pling. Among the relevant operators is detuning v from 
the Feshbach resonance in the absence of the lattice, and 
the co-efficient that contains it is fixed by relating the 
scattering matrix of this theory to the scattering length 
a = -ljv for V(r) = 0. 

Our primary result is that the above model has a uni- 
versal phase diagram as a function of fJ,/E z , E T /V, and 
uajj, where E x = tt 2 h 2 / (Ama^) is the molecular recoil 
energy. The phase boundaries are shown in Fig. Q] as a 
function of the first two parameters for different values 
of va^ ■ 

We will explore the superfluid phase boundary by fo- 
cusing on the superfluid order parameter. We integrate 
out the fermion fields ipi& and obtain the effective action 
of the Cooper pair field which can be expressed using 
Feynman diagrams: 



= N— I drdW$ 

47T 




(3) 



Each fermion loop (straight lines) contributes a factor of 
N, so that N appears as an overall factor in the effec- 
tive action. Every external wavy line represents the pair 
field $ or $t, depending on the direction of arrow with 
respect to vertex, and every vertex is associated with a 
point in space and time. In the present problem, the op- 
tical lattice potential V(r) breaks translation symmetry, 
so that momentum is not a good quantum number. In- 
stead, according to the Bloch's theorem, quantum num- 
bers are band index n and crystal momentum fc which 
takes values within the first Brillouin zone (BZ). Assum- 
ing a Bravais optical lattice with reciprocal vectors G, 
the bare fermion propagator is given by: 



G (l,2) 



BZ 



x e 



i(fe(r 2 — ri)— o>(t 2 — n)) i(G 2 r 2 — Giri) 



(4) 



where the Fourier components of the Bloch wavefunctions 
ipn,k-,G and energies t n ,h are obtained from the Fourier 



transformed Schrodinger equation: 
'(k + G) 2 
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G' 
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V G - G > 



4>n,k;G' 



(5) 



Solving this set of equations numerically requires impos- 
ing a cut-off energy A that limits the reciprocal vectors 
kept in calculations: G 2 < 2mA. For our V(r) above, 
the only non-zero Fourier components are Vq = 3V/2, 
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The presence of superfluidity can be described by a su- 
perfluid order parameter, which due to the optical lattice 
may have Fourier components $g at various reciprocal 
lattice vectors G. Fluctuations S&(r) of the Cooper pair 
field are added to the order parameter, and the total bo- 
son field $(r) = ^ G $Ge lGr + 6$(r) is represented by 
the wavy lines in J3J . For large N the action S c g oc N be- 
comes large, so that fluctuations are suppressed; the 
mean- field theory becomes exact in the limit N — > oo. 
Integrating out <5$ gives rise to corrections of the order 
1/N to the mean-field results. This follows from dia- 
grammatic perturbation theory performed in the effec- 
tive action, where the first diagram in ([3]), together with 
the detuning term, defines a 'bare' propagator of the $ 
fields, while all other diagrams define new vertices of the 
$ fields. Even though the new vertices are proportional 
to N, the 'bare' propagator is proportional to 1/N and 
hence yields perturbative expansions of thermodynamic 
functions in powers of 1/N. 

As in the bosonic case, we expect a second order 
superfluid-insulator transition. Near such a transition, 
the action terms quadratic in order parameter determine 
the state, in analogy to a simple uniform $ 4 Landau- 
Ginzburg theory of bosons. Neglecting the fluctuations 
5Q(r) in the N — > oo limit, the free energy density T is 
just the effective saddle-point action divided by volume 
V and $ = 1/T: 



N 



•err 



N/3V 



( 2 ) St 



G'^G 



+ G($ 4 ). (6) 



G.G' 



(2) 

The quadratic couplings K G ' G , can be represented as 
a matrix whose components are indexed by the recip- 
rocal lattice vectors G. If this matrix has only posi- 
tive eigenvalues, the free energy ([6]) is minimized by 
(VG) $g = indicating an insulating (T = 0) or normal 
(T > 0) phase. Otherwise, the minimum is obtained at 
(3G) $g 7^ 0, and the established phase is superfluid. 
The role of 0(<E> 4 ) terms near a second order phase tran- 
sition is only to stabilize the theory. 

A naive derivation of K G G , from the first Feynman di- 
agram in (|3|) produces an ultra-violet divergent expres- 
sion. This divergence stems from the naive continuum 
form of the bare field theory, and must be renormalized 
away by absorbing it into finite physically measurable 
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renormalized quantities. One step toward this goal has 
already been taken by absorbing any bare molecule mass 
(a part of the bare molecule dispersion omitted from ©) 
into detuning v 1 which is measurable and fixed in the ef- 
fective field theory simply by the properties of scatter- 
ing matrix. The second step is to remove the remaining 
unphysical divergent part by dimensional regularization, 
which at large momenta k + G is carried out just like in 
a system without the optical lattice @. The regularized 
expression is: 



0.075 



G,G' = i^G,G' 
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E E [^-f(^uk)-f(e n2 ,-k)]x (7) 
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where f(x) = {l+e x l T )~ 1 is the Fermi-Dirac distribution 
function. 

In the following, we numerically compute the matrix 
f or reciprocal lattice vectors \G\ < V2toA', where 
A' < A, in the limit N — > oo. In practice, it is sufficient 
to choose a very small A' (\G\ < 2na^ 1 ) 1 as long as A 3> 
V . By mapping dependence of the smallest eigenvalue 
I^ 2 ) of the matrix on chemical potential fi, lattice 
amplitude V and detuning v, we find the second order 
phase boundary between superfluid and insulating phases 
by the condition 



r« = o. 



(8) 



In FigfTJwe show a contour plot of the phase boundary as 
a function of chemical potential and inverse lattice am- 
plitude at T = and N — oo. This choice of plot axes 
was made in order to obtain resemblance to the well- 
known phase diagram 16] of a superfluid to Mott in- 
sulator transition; while chemical potential directly con- 
trols density, inverse lattice amplitude is correlated with 
hopping strength t of an effective tight-binding model. 
Indeed, insulating regions are dome-shaped and corre- 
spond to integer fillings of optical lattice sites with bosons 
(Cooper pairs) . The character of an insulator depends on 
detuning v from the Feshbach resonance. On the BCS 
side of the resonance, v > 0, the insulating domes be- 
come larger, and converge toward the fermion band-gap 
boundaries as v is increased. The system then behaves as 
a typical band-insulator (weakly paired fermions are in- 
sulating due to a filled band). On the other hand, in the 
BEC limit v < the domes become 'smaller', resembling 
a Mott insulator. When molecules are tightly bound, 
filling up a fermion band is not sufficient to destroy su- 
perfluidity, but repulsion between molecules needs to step 
in. Note, however, that true Mott insulating phases, with 
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FIG. 1: A contour plot of the superfluid-insulator phase 
boundaries at T — and N = oo. Contours, dependent on 
detuning v from the Feshbach resonance, separate superfluid 
regions surrounding the fermion bands (black) from insulat- 
ing regions in the band structure gaps. The dashed yellow 
contours correspond to the transition at the resonance v — 0, 
while the spacing between contours is Au = a^ 1 , where ol is 
the optical lattice spacing; contours move upward and toward 
the band edges as v grows. Insulators are labeled by the 
closest even integer to the average filling of the optical lat- 
tice by atoms in the grand-canonical ensemble. The reference 
energy scale is molecule recoil energy E T — it 2 ft 2 / (Ama^) . 
Accuracy is smaller than that in Fig[2] because calculating 
at many points in reasonable time required a small cut-off 

A = 3 x (27TOL 1 ). 



an arbitrary integer lattice filling by molecules, cannot be 
found without including 1/7Y corrections. 

It is also important to consider the phase diagram as 
a function of particle density, rather than chemical po- 
tential. The Mott insulating lobes in Fig. [1] all have a 
density of an even integer number, np, of fermions per 
unit cell. We therefore fix the density at n-p and study 
the transition from the insulator to the superfluid. By 
generalizing the argument made for the bosonic case, this 
T = transition occurs at the point where we satisfy the 
condition ([8]) along with 



dp, 



(9) 



The two conditions (jHJ) and © determine an isolated 
point in the phase diagram of Fig. [1] for each np, and the 
location of these points then immediately yields Eq. |T]). 
In Fig. [5] we plot the universal function Fn-p^va^) for 
n-p = 2 and N = oo (calculations were performed with 
large cut-off A < 10 x ^Tra^ 1 )). 

The obtained values of V c = F nF (ai J i')E r are larger 
than those reported in the experiment, V c w 6E Y [lOj. 
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FIG. 2: Universal function F nF (iyah) in {TJ at the superfluid- 
insulator transition for the first fermion band completely 
filled, n F = 2 (T = 0, N = oo). 




FIG. 3: A contour plot of the superfluid-normal phase bound- 
ary atT^ 1.8-Er and N = oo. This temperature is T w 0.6i5f, 
where _Ef = fi, 2 (67raL) /(2m) is the Fermi energy of a free 
fermion gas at the same density as the fermion gas in the lat- 
tice with two particles per site. The black and yellow contours 
are contours of constant detuning, as in Fig[T] The green lines 
are normal-phase constant density lines, starting at tif — 1 at 
the left-bottom, displaced by AriF = 1. The superconductor- 
normal phase transition can occur at any density or detuning, 
at the intersection of the corresponding contour lines. Note 
that due to a small cut-off used in plotting this diagram, one 
should not directly compare it with Fig. [5] - the purpose is 
only to illustrate the shape of phase boundaries. 

However, unlike our computation so far, the experiment 
is at a non-zero temperature (not known very accurately 
in the presence of the optical lattice), and thermal fluc- 
tuations should decrease the value of V c . We extended 
our results to a range of T > as shown in Fig. [3] The 
character of the transition at which superconductivity is 
lost changes qualitatively at T > 0. At T = 0, the in- 
sulator must have a density rip = even integer, and so 



a superconductor-insulator transition can only occur at 
such values of rip; this was rationale behind the addi- 
tional constraint in Eq. ©. However, at the T > 0, 
the transition is more properly a superconductor-normal 
transition, and the normal state can occur at any den- 
sity. In Fig [3] we show a contour plot of the boundaries 
between supcrfiuid and normal phases at finite temper- 
atures (N — oo). As temperature is increased, the non- 
superfluid domes gradually expand. The contours cor- 
responding to larger values of v (BCS-limit) are more 
affected by thermal fluctuations than those correspond- 
ing to smaller values of v (BEC-limit). Since at T > 
the normal regions can occur when the chemical potential 
is not in a band gap, the normal regions corresponding 
to different average lattice fillings can merge when the 
independent effect of fluctuations is large enough; when 
this happens at a particular v, the appropriate contour 
stretches all the way from the bottom to the top of the 
diagram, instead of being dome-shaped. We expect that 
including molecule fluctuations would increase these ef- 
fects even further. Now, if an experiment is performed at 
a fixed density (green lines in Fig.[3J, the phase transition 
to a normal phase at unitarity (dashed line) can occur at 
a smaller lattice depth V c than at T — 0. The reduc- 
tion of V c can be particularly dramatic if the transition 
is observed in a region where the T = non-supcrfluid 
regions have merged due to thermal fluctuations. This 
can even occur at temperatures small compared to the 
Fermi energy. We suspect this effect is the primary rea- 
son for the discrepancy between the T > experiments 
and our T = results. 
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